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Abstract 

The new particle X recently discovered by the ATLAS and CMS Collaborations in 
searches for the Higgs boson has been observed to decay into 77, ZZ* and WW*, but 
its spin and parity J^, remain a mystery with = 0"*" and 2+ being open possibilities. We 
use PYTHIA and Delphes to simulate an analysis of the angular distribution of gg ^ X ^ 
decays in a full 2012 data set, including realistic background levels. We show that this an- 
gular distribution should provide strong discrimination between the possibilities of spin zero 
and spin two with graviton-like couplings: ~ 3(j if a conservative symmetric interpretation 
of the log-hkehhood ratio (LLR) test statistic is used, and ~ 6(j if a less conservative asym- 
metric interpretation is used. The WW and ZZ couplings of the Standard Model Higgs 
boson and of a 2+ particle with graviton-like couplings are both expected to exhibit custo- 
dial symmetry. We simulate the present ATLAS and CMS search strategies for X — > l^W* 
using PYTHIA and Delphes, and show that their efficiencies in the case of a spin-two particle 
with graviton-hke couplings are a factor ~ 1.9 smaller than in the spin-zero case. On the 
other hand, the ratio of X2+ WW* and ZZ* branching ratios is larger than that in the 
0+ case by a factor ~ 1.3. We find that the current ATLAS and CMS results for X WW* 
and X — > ZZ* decays are compatible with custodial symmetry under both the spin-zero and 
-two hypotheses, and that the data expected to become available during 2012 are unlikely 
to discriminate significantly between these possibihties. 
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1 Introduction and Summary 



A new particle X with mass ~ 125 to 126 GeV has been discovered by the ATLAS [T] and 
CMS f2] Collaborations during their searches for the Higgs boson of the Standard Model. At 
first sight, the new particle X is observed to have similar characteristics to the long-sought 
Higgs particle H: it is a boson that does not have spin one, and hence it is different in 
character from gauge bosons. But is it really a (the) Higgs boson? 

Answering this question will require a number of consistency checks. For example, are 
the X couplings to other particles (at least approximately) proportional to their masses? 
A first step towards answering this question was taken in [s] (see also jij), where it was 
shown that the data published in the ATLAS and CMS discovery papers (T||2] and other 
public documents, in combination with results from the TeVatron experiments CDF and 
DO [5], are inconsistent with mass-independent couplings to the t, and 6, and in fact 

highly consistent with linear scaling ~ M/v, where v ~ 246 GeV is the expected electroweak 
symmetry-breaking scale. This type of consistency check will be improved significantly with 
upcoming data. 

However, an even more basic question is whether the recently-discovered Higgs candidate 
has the spin-parity = O"*" expected for the scalar Higgs boson in the Standard Model, 
and several ways to test this have been proposed. For example, some of us have recently 
pointed out |6j that the XV invariant mass distributions in associated production of the 
X particle together with a vector boson V = Z or W are in theory very different for the 

assignments O"*", 0~ and 2"*" for the X particle (where we assume graviton-like couplings 
in the last case). We have also shown that these differences in the VX mass distributions 
are maintained in simulations with realistic detector cuts applied, and hence may be used 
to obtain indications on the of the X particle, if the experimental backgrounds can be 
suppressed sufficiently. At the time of writing, the TeVatron experiments CDF and DO have 
reported evidence for X production in association with vector bosons V, followed by X — )■ 66 
decay, at a rate compatible with the Standard Model [5j, but have not yet reported the VX 
mass distributions. The most important backgrounds are expected to have small invariant 
masses. (Non-)observation of a similar VX signal in ATLAS or CMS at the Standard Model 
level with an invariant mass distribution corresponding to the 2+ prediction would provide 
evidence for (against) the 2+ hypothesis. 

Here we first study the potential discriminating power of the angular distribution of 
(yfgf — )■ X — )■ 77 events, which we simulate using PYTHIA and the Delphes, including realistic 
background levels. We find that the data already available may be able to offer some dis- 
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crimination, and that the data hkely to be available to the ATLAS and CMS Collaborations 
by the end of 2012 should be able to distinguish between the spin- zero and graviton-like 
spin-two hypotheses at about the 3-0" level if a conservative symmetric interpretation of the 
LLR test statistic is used, or over 6 cr if a less conservative asymmetric interpretation is used. 

There have been many suggestions to exploit correlations among decay products of the 
'Higgs' to identify its spin and parity [?]. Some of these correlations have been incorporated 
in the Higgs search strategies adopted by the ATLAS, CMS, CDF and DO collaborations. 
These are based on the assumption that its = 0"*" and have, in general, lower efficiencies 
for detecting a particle with different J^. For example, the ATLAS and CMS searches for 
H — )■ WW* — )■ decay make use of the kinematic correlations expected for a scalar 

particle, and the search by CMS for H — t- ZZ* — )■ Ai^ also exploits the correlations expected 
in decays of a spin-zero particle. 

We study here whether the logic can be inverted to argue that the results of these searches 
already favour the spin-zero hypothesis for the X particle over the spin-two hypothesis. We 
find that the current ATLAS and CMS measurements favour custodial symmetry for the 
XWW and XZZ both if the X particle has = O"*" and if it has = 2^ and graviton-like 
couplings. This result is based on simulations of the ATLAS and CMS H — )■ WW* searches 
using PYTHIA and Delphes, which indicate that their efficiencies are a factor ~ L9 lower 
under the spin-two hypothesis. On the other hand, this effect is partially offset by the ratio 
of the X — )■ WW* and ZZ* branching ratios that is larger by a factor ~ 1.3 in the 
case. Extrapolation to the full expected 2012 data set suggests that it will not be able to 
discriminate significantly between the O"*" and graviton-like 2"*" hypotheses using the ratio of 
X WW* and ZZ* decays. 

2 Possible spin-parity assignments 

As is well known, the fact that the new particle has been observed to decay into a pair of 
on-shell photons implies that it cannot have spin one [l|^- Accordingly, here we consider 
the spin-zero and spin-two options. In the spin-zero case, we consider the pseudoscalar 
possibility = 0~ as an alternative to the assignment O"*" expected for the Standard Model 
Higgs boson, and we consider spin-two models with graviton-like couplings. 
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2.1 Pseudoscalar couplings 

In this case, the following are the couplings to two vector bosons and a fermion-antifermion 
pair, respectively: 

V;K: ^,up.v'q% If-- 75- (1) 

where the first term corresponds to an e^yp^F^'^ F'^" X interaction term, with and the 
sum and difference of the four-momenta of the two vector bosons, respectively. In this case, 
the forms of the vertices are unique, though the normalizations are arbitrary. We assume a 
custodial symmetry so that the pseudoscalar couplings to and Z are equal, but make 
no assumptions about the ratio of the pseudoscalar couplings to the photon and gluon. 

2.2 Tensor couplings 

Several forms are possible for the couplings of a spin-two particle to two vector bosons. It 
was shown in js] that Lorentz and Standard Model gauge symmetries forbid couplings of 
a massive spin-2 particle to two Standard Model particles through dimension-four terms in 
the Lagrangian and, assuming the flavour and CP symmetries of the Standard Model, it 
should couple flavour-diagonally via dimension-flve terms that take the same forms as their 
energy-momentum tensors, namely 

C.nt = -^G^^Tl,, (2) 

where the T^^ are the four-dimensional stress tensors of the Standard Model species i = 
b,f,V,..., where V denotes a generic gauge boson. In scenarios with extra dimensions, 
Meff is the effective Planck mass suppressing the interactions {M^jf — O(TeV)), whereas in 
composite models it denotes a scale related to conflnement. 
In the speciflc cases of vector boson flelds V^, one has: 

= -F^Fpy + ^ z/) - mlV.Vy , (3) 

where F^^y is the fleld strength tensor for and the vector boson mass term oc my would 
be absent for the photon and gluon. In the case of a fermion /, one has 

f i - - 
T^u ^ 2 ^^^^''^ + (/i ^ z/) + h.c. -nifi)i)gpy . (4) 

Since the couplings of a composite spin-two particle and a massive graviton Standard 
Model particles would both take the forms ([2]), the model dependence would appear in the 
coefficients q. In the case of a resonance in a strong sector, these coefficients would reflect 
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the underlying dynamics and the quantum numbers of the constituent fields. For example, 
if the constituents do not carry color, the coupling to gluons, Cg, would be zero, whereas the 
coupling to photons, c^, would refiect their electric charges. 

A more specific scenario is that of massive gravitons in warped extra dimensions, with 
the Standard Model particles residing in the bulk j9|. In five-dimensional scenarios with a 
factorizable metric 

ds'^ = w'^{z){rif,^dx''dx'' - dz"^) , (5) 

one may consider various possibilities. For example, whereas in a fiat extra dimension 
w{z) = 1, in an AdS extra dimension w{z) = zjjy/z. If a Standard Model field lives on a 
brane located at z^ G {zuviZir), one has 

c-^. (6) 

In fiat extra dimensions there would be no parametric suppression, w = 1 and the couplings 
(|6| would be universal However, in warped extra dimensions with w{zir) <^ w{zuv) and 
for a field living on the UV brane c = w{zjji) /w{zuy) ~ 1 TeV/Mp [Toj. On the other hand, 
couplings to massless gauge bosons would be suppressed by the effective volume of the extra 



dimension, namely 11 



c ~ 1/ / w{z)dz . (7) 

J zuv 

In the AdS case, this suppression would be by log^zm/zuv) — log(Mp/TeV) ~ 30, whereas 
in fiat space the suppression would be by the entire volume of the extra dimension. In either 
case, the couplings for gluons and photons would be equal: Cg = c^, assuming that localized 
kinetic terms are not unnaturally big. 

We note that, in either case, custodial symmetry is ensured for the spin-two couplings 
to the massive vector bosons of the Standard Model. The reason is gauge invariance: the 
graviton couples to the gauge eigenstates W"" universally, which implies that cw = cz as 
g' — )■ 0. Once electroweak symmetry breaking occurs, the graviton would feel the effect 
through couplings like Gfj^uD^T,^ D'^'E, where S is a physical or spurion field, and (S) = v. 
As S should respect an approximate custodial symmetry (as indicated by the small value 
of the T and Ap parameters (l2]), the spin-two particle inherits the approximate custodial 
symmetry. 

If we assume that electroweak symmetry is broken by boundary conditions on the IR 
brane, we expect that the support of the W and Z wave functions will be suppressed in the 



"'^We note that in flat extra dimensions one typically finds conservation of a Kaluza-Klein parity that 
would forbid the coupling of the massive graviton to two light Standard Model particles. 
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neighbourhood of this brane, so that cwt,Zt < c-y,g, where the Vt are the transverse compo- 
nents. On the other hand, the longitudinal W and Z are localized near the IR brane 13 , as 
the wave functions of massive fermions such as the b and t, so that cv^^^^^, Cb^t > c^,g- and 
the wave functions of light fermions such as the u and d are expected to to be concentrated 
closer to the UV brane, so that Cu^d ^ c^,g- 

To summarize, in warped extra dimensions one expects the following qualitative be- 
haviour 

Warped AdS : Cb, Ct > cw — cz ^ = Cg Cu, Cd ■ (8) 

with, e.g., Cb — 30c^^g in the case of a Randall- Sundrum model with the third generation 
located near the UV brane [9 , 14 . In Section 4 below, we focus on the question whether the 



custodial symmetry relation cw — cz is compatible with the present experimental data. 

3 Angular distributions in Xq 2 — > 77 decay 

In the case of a spin-0 particle, Xq, its decay products are isotropically distributed over a 
two-dimensional sphere, so one expects a flat distribution as a function of cos^*, where 6* 
is the angle of the photon relative to the beam axis in the X rest frame. On the other 
hand, the 77 angular distribution will in general be non-isotropic in the case of a spin-2 
particle, X2. Assuming that the gluon-gluon fusion process dominates X2 production, and 
graviton-like couplings of X2 to both gluons and photons as discussed in Section 2, the 77 
angular distribution in the X2 centre of mass in gluon-gluon collisions was calculated in |15j 
(see also [16]), and is given by 

l^i^cos^r.lcos^r , (9) 

which can in principle be distinguished from the isotropic distribution expected in the O"*" 
and 0~ cases. 

In order to see whether these theoretical distributions are distinguishable in practice, we 
simulated samples of spin-0 and spin-2 production accompanied by 0, 1 or 2 jets, followed 
by decay to two photons: 

pp^Xo,2(^77) + (0,l,2) jets (10) 



using MadgraphS 17 , which implements graviton-like couplings in the 2"*" case. Our simu- 



lation includes gluon fusion, vector boson fusion and production in association with the top 
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and vector bosons, in the same proportions as in the case of a conventional Higgs boson |^ 
The events are then matched using the MLM scheme in PYTHIAv6 . 4 18 , and passed through 
the Delphes [l9) simulation code. 



Figure [T] displays the cos 9 (left) and cos 9* (right) distributions after implementing the 
baseline cuts > 20 GeV and \ri\^/ < 2.5. We see that the theoretical difference between 
the scalar Higgs and graviton-like decay distributions in the rest frame of Xo,2 survives 
these basic cuts. 
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Figure 1: Simulation of the Xo_2 — ^ 77 polar angle distributions obtained using MadGraphS, 
PYTHIAv6.4 and Delphes, after implementing the baseline cuts > 20 GeV and \ri^\ < 
2.5. The left panel show distributions in the laboratory frame, and the right panel show 
distributions in the diphoton centre- of-mass frame. 

We have studied whether the higher-level selection cuts could affect the distributions and 
the discriminating power between the spin-zero and -two hypotheses. As shown in Fig. [2| we 
find that the distinction is quite stable under changes in the photon momentum cuts, e.g., 
to p'^'^ > 40, 25 GeV, or in the pxt cut separating the glue-glue and vector-boson-fusion- 
enhanced processes 

3.1 Toys and Statistical Procedure 



As a first step in our analysis, we use a simple angular asymmetry variable as in 20 



to 



quantify the separation significance between spin O"*" and 2+ as a function of the number of 



^To the extent that these other production mechanisms are suppressed relative to gg — > X, their inclu- 
sion or omission is not important. We have checked that their inclusion in our simulation docs not affect 



significantly the angular distributions from gg ^ X alone. 

•^We define prt = Ipt' ^ where the thrust vector is defined as t = {p^ — Pt)/\Pi 
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Figure 2: The effects of different selection cuts in p^'^ and pxt on the angular distributions 
for the spin-zero case (left panel) and the graviton-like spin-two case (right panel). 



signal events in Monte-Carlo (MC) simulations, starting with the idealized case of only signal 
events, then showing how the asymmetry can be extracted in the presence of background 
and describing how this can be simulated with toy MCs. The result is presented for different 
signal to background (S/B) ratios representing different event sub-categories. We then repeat 
the analysis using the complementary log-likelihood ratio (LLR) test statistic, finding results 
that are somewhat more sensitive B 

A reference sample of 10000 spin-zero signal diphoton events from the process pp h ^ 



77 was generated using MadGraphS vl.4.8.3 17 and passed through an ATLAS detector 



simulation based on Delphes [i9j. After transforming the diphoton system to its centre-of- 
mass frame, baseline Pt cuts of 40 and 25 GeV were applied on the leading and sub-leading 
photons, respectively. As shown above, the angular distribution does not vary appreciably 
before and after the detector simulation and cuts. The angular distribution of this reference 
sample was reweighted to obtain a spin-two reference sample. These reference histograms 
were then sampled repeatedly to provide toy histograms with numbers of signal events that 
could be expected realistically. We have checked that this procedure gives results similar to 
generating each toy individually. 

For each toy, we first quantify the shape of the distribution in cosO* by an asymmetry 
variable, defined as 

A^centre Apsides 



A 



11 



A^centre ~l~ Apsides 

where A'centre is the number of events lying within the range —0.5 < cos9* < 0.5 and Apsides 



For a discussion of the relative merits of the LLR and the asymmetry variable, see 20 . 
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is the number of events outside this range. Populating a histogram of the asymmetry value 
for each toy gives a distribution around different means for the spin-zero and -two toys, as 
illustrated in Fig. [3] for 10000 toys of 160 signal events. 




Figure 3: Distribution of the signal angular asymmetry variable A (11) for 10000 toys of 160 
signal events each, with a superimposed Gaussian fit in red. The histogram for the spin-zero 
toys is unshaded, and that for the spin-two toys is shaded blue. These plots do not take 
backgrounds into account. 



Using the asymmetry (11) as our test statistic, with a distribution that is fit well by a 



Gaussian (as seen in Fig. |3]) to obtain a normalized probability distribution function pdf(A), 
we quantify the separation significance using two different methods denoted as 'asymmetric' 
and 'symmetric'. 

In the asymmetric method one is biased towards verifying hypothesis Si and excluding 
hypothesis 5*2. Thus the value of the asymmetry that is expected to be measured by an 
experiment is taken to be the mean of the distribution for 5*1, namely A^^*^*^ = A^'^^^, and the 
extent to which we can exclude hypothesis S2 is the area (3 under the tail of pdf52(A): 



jj^ pdf^^(A)rfA 
rooPdf5.(A)rfA 



For A^'^^^ > A^^'^^ the integral is taken over the other side of the distribution, so as to obtain 
instead the area under the right tail of the distribution of 82- The quantity P is also known 
as the 'Type II' error, namely the probability of wrongly accepting hypothesis A in the case 
that actually 5*2 is true. 
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In the symmetric method the two hypotheses are treated equally. One defines a Acut-ofr 
for which a = /3, namely the area a under the right tail of 5*1 is equal to the area /3 under the 
left tail of 5*2 (in the example where A^^^^ < A™^'^"). Thus, whatever the value of Aobs found 
by an experiment, if it lies to the left of Acut-os hypothesis Si is accepted and 5*2 rejected, 
and vice versa if Aobs > Acut-ofr- The significance is given hj a = (3. 

Both approaches are justified in that there is strong motivation for prioritizing the spin- 
zero hypothesis, and thus quoting the asymmetric significance (see also [20j). On the other 
hand, the symmetric approach is more objective and conservative (see [2l]). In the following 
we quote results obtained using both methods. 

The significance a is translated into n standard deviations by finding the equivalent area 
under a standard Gaussian distribution B 

a = I e'^^dx . (12) 



For example, a = 0.05 corresponds to n = 1.64, and the discovery standard of n 
corresponds to a = 2.87 x 10"''. 



Asymmetry Variable, Signal Only 
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asymmetric c 



Figure 4: Numbers of standard deviations a of separation (12) for two different methods 
of interpreting the statistical significance, as functions of the numbers of signal events with 
no background taken into account. The blue (red) dots are obtained using the asymmetric 
(symmetric) method, respectively. The dotted lines on the left and right are the expected 
yields in high and low S/B categories, respectively, with 30 fb~^ of integrated luminosity. 

""This is the one-sided definition most commonly used in the hterature, as opposed to the two-sided 
convention sometimes seen, which generally yields a higher number of standard deviations for the same 
p- value. 
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In Fig. |4] we show the separation significance defined in this way for the symmetric 
and asymmetric interpretations as functions of the numbers of signal events, neglecting 
backgrounds. The vertical dotted lines indicate the expected number of signal events in 
two categories that can be reached with 30 fb~^ of 8 TeV data collected by the end of the 
year. The expected yields in the two diphoton categories are obtained by combining the 
CMS BDT categories 0,1 and 2,3 from Table 2 of [2]. These correspond to high and low 
signal-to-background (S/B) ratios of approximately 0.42 and 0.19, respectively. 

Note that we have not included backgrounds in Fig. |4} This figure should therefore be 
taken as an idealized limit of what could be achieved assuming a perfect separation of the 
signal from the background. In the next Section we give an estimate how the background 
would affect this. 

3.2 Background Simulation 

A reference sample of 20K QCD continuum background pp — )■ 77 events with the parton- 
level invariant masses of the diphoton pairs between 124 and 126 GeV was generated. The 
same detector simulation and cuts as above were applied. This sample was then used to give 
a number of toy background events corresponding to the desired S/B ratio. 

For each toy the background and signal events are added together to give a total sample 
representing the available experimental information. All that can be measured is the total 
asymmetry in the signal region, Atot- However the signal asymmetry can be extracted, since 

I A^centre Apsides 

^ iw:+A^ 

( Af^ -L ]\[^ ] — ( Af^ -I- /V^ 1 
V centre ' centre/ V sides ' sides/ 

centre ' ^ centre i sides ' sides 

= fA, + {l-f)A, , (13) 

where we assume that / = Ns/{Nb + Ng) = (Atot — A'fe)/Atot is known, and the asymmetry 
of the background, Af,, can in principle be measured with high accuracy in the invariant 
diphoton mass sideband regions outside the signal region. 

The error in the experimental determination of / is the main limiting factor in recon- 
structing Ag. To simulate the effect of this, for each toy we randomly sample A";, from a 
Gaussian centered around the true value of A";,, with a one-a width of a/A^, the statistical 
error. We assume that Atot is measured much more accurately, so that its error can be 
neglectecPj We calculate A^ using the background reference sample, so as to simulate the 



^Adding a 1% error in A^tot does not affect the results. 
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measurement from the sidebands that is assumed to have much higher statistics. 




Figure 5: Distribution of the extracted signal angular asymmetry variable A (11) for 10000 
toys each. The histogram for the spin-zero toys is unshaded, and that for the spin-two toys 
is shaded blue. These simulations include the background with values of S/B = 0.42 (left) 
and 0.19 (right). 



As mentioned above, the benchmark luminosity that the experiments hope to attain in 
2012 is 30 fb~^ at 8 TeV. The CMS diphoton search separates events into categories that can 
be approximated by two samples with high and low S/B = 0.42 and 0.19, respectively, with 
~160 and 420 expected signal events respectively for 30 fb^^. The resulting distribution of 
As extracted from this simulated measurement of Aot? ^b, and / per toy is illustrated in Fig. 
5] for 160 signal events in the high S/B category. Fig. |6] shows the separation significance as 
a function of luminosity for the two categories. 

We see from Fig. [6] that these simulations translate into a significance of 3 to 3.5 a in 
the asymmetric interpretation for each category. Ideally, a more detailed simulation should 
be done in which the events for each category are output from the BDT that sorts them, as 
this may affect the angular distribution more significantly than the simple cuts used here. 
As a basic check, we have verified that placing \ri < 1.4| cuts on both photons to simulate 
barrel-barrel events does not alter our results substantially. 
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Figure 6: Numbers of standard deviations a of separation (12) for two different methods of 
interpreting the statistical significance of the measurement of the angular asymmetry A, as 
functions of luminosity with S/B = 0.42 (left) and 0.19 (right). The blue (red) dots are 
obtained using the asymmetric (symmetric) method, respectively. Dotted lines indicate the 
luminosity expected by the end of 2012. 

3.3 Log-Likelihood Ratio Test Statistic 

An alternative test statistic is the log-likelihood ratio (LLR), for which the likelihood Cs for 
spin hypothesis in a single toy pseudo-experiment is defined as 

events 



^5 = n p^f^ia^i) 



The probability density function in cos 9* of the signal is extracted from the probability to 
lie in a bin of the high-statistics MC from which the toys are sampled. After calculating the 
likelihoods for both hypothesis Si and 5*2 we obtain the LLR for that toy by taking —2 In 



'S2 



Thus, if we generate a set of toys for spin hypothesis 5*1 (52), the LLR distribution will be 
centred around a negative (positive) mean. We may the quantify the separation between 
these two distributions as described previously. 

For a pure signal with no backgrounds, the LLR distributions for spin O"*" and 2"^ with 
160 signal events are shown in Fig. [7| Also plotted in Fig. [7] is the separation significance 
as a function of the number of signal events, though we emphasize that this is in idealized 
limit in which the signal events can be perfectly extracted from the background. 

In order to include the effects of the background, we combine the background and signal 
MC to form a total MC, then simulate the extraction of the signal by subtracting statistically 
the number of background events expected in each bin from the total number of events in 
that bin. The number of background events per bin is smeared using a Gaussian centred on 
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Figure 7: An example distribution of the test statistic for 160 and 430 signal events, above, 
and the separation significance obtained using the LLR test statistic as a function of the 
number of signal events, bottom. Note that no backgrounds are included here. 



the true value, with a one-a width given by the statistical error y N^^. If the randomly 
smeared number of background events exceeds the total number of events in that bin the 
corresponding bin of the measured signal histogram is set to zero (since it cannot be negative). 

Distributions in LLR and separation significance plots similar to those in Fig. [7] are 
shown in Fig. [8} with the backgrounds now taken into account, for both high and low values 
oi S/ B = 0.42 and 0.19, respectively. We see that the high (low) S/ B category for 30 fb~^ of 
8 TeV data, corresponding to 160 (420) signal events, yields a separation significance around 
3.5 (4.2) cr using the asymmetric method. A combination of the high and low categories can 
achieve over 6-cr separation using the asymmetric method, as is seen in Fig. |9j and ~ 3 a 
using the symmetric method. 
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Figure 8: Separation significance using the LLR test statistic as a function of the lumi- 
nosity (lower panel) and an example distribution of the test statistic for the 160 and 420 
signal events expected in the high- and low-S/B categories at 30 fb~^ (upper panels). The 
backgrounds are included, with S/B = 0.42 (left) and 0.19 (right), respectively. 

4 Angular, mass and rriT Distributions in Xq 2 — > WW* — > 
Decay 

As discussed in ^6], the lepton momentum distributions and correlations are very different 
for the Xq^2 hypotheses. In the spin-zero case, the spins of the and W^* must be antipar- 
allel, implying that the charged leptons £^ produced in their decays appear preferentially in 
the same hemisphere. On the other hand, in the spin-two case, the spins of the and 
W^* must be parallel, implying that their daughter appear preferentially in the opposite 
hemisphere. As pointed out in |l6] , these differences in the decay kinematics imply that the 
dilepton invariant mass m« is generally smaller in Xq decay than in X2 decay, as is the dif- 
ference (pee between the azimuthal angles. A corollary is that the net transverse momenta 
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Figure 9: Separation significance using the LLR test statistic as a function of luminosity 
for the combination of high- and low-S/B categories. The dotted line indicates the expected 
reach by the end of 2012. 



of the pair, p^'^^, is generally larger in Xq decay than in X2 decay. We now address 
the question whether and to what extent these differences survive the event selections and 
cuts made by ATLAS and CMS. 

Regarding the simulation details, we created new models in Feynrules [22j, including the 
pseudoscalar and graviton-like spin-two couplings described in Section 2. We then interfaced 
with MadGraphS 17 using the UFO model format 23 . We incorporate hadronization and 



showering effects using PYTHIA 18 , and detector effects using Delphes 19 . 



4.1 Simulation of ATLAS and CMS event selections 

We start by implementing the baseline cuts of the ATLAS and CMS analyses: two isolated 
leptons (= e, /i) of pt > 15 GeV and |?7| < 2.5. As seen in Fig. 10, the distributions in the 
dilepton invariant mass m« (left) and in the relative azimuthal angle (pu (right) are very 
different in the baseline Xq and X2 simulations. These differences reflect the kinematical 
effects noted earlier. In particular, ma is generally smaller in the Xq case than in the X2 
case, as seen in Fig. 7 of 16 , as is reflecting the angular distribution shown in Fig. 6 
of [16 . 
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Figure 10: Results from simulations oftheXQ+ 2+ WW* — )■ l^vil'vi signals, using PYTElk 
and Delphes and the baseline cuts described in the text. The left (right) panel displays the 
(4>u) distributions. 

Our next step is to simulate the ATLAS search for X — )■ WW* — events. 
This is based on a selection of events with two opposite-sign, unlike-flavour leptons with 
p^'^^ > 25, 15 GeV in the central region, and invariant mass mu G [10,80] GeV0 The 
events are then separated into categories with 0, 1 and 2 anti-Zc^ jets (defined by cones 
R = 0.4) and pt > 25, 30 GeV in the central and forward regions, respectively. In the 0- and 
1-jet samples used here, the dilepton invariant mass upper bound is tightened to 50 GeV. 
The following set of cuts is then applied to the 0-jet sample: 

E^p^i > 25 GeV, > 30 GeV, and |A0«| < 1.8 , (14) 

where EJj?'' = pj- sin A0mj„ with A^^m = min(A0, vr/2) and A0 the minimum angle between 
the missing-energy vector and the leading lepton, the subleading lepton or any jet with pt > 
25 GeV. In the 1- and 2-jet case, there is an extra cut 

\ = \f^+p)^ + I < 30 GeV ( 15) 

as well as a 6-tag veto. In the 2-jet case we also implemented the vector-boson-fusion cuts 

rrijj > 500 GeV and \ Ayjj\ > 3.8 . (16) 

Finally, a cut = a/ {E^ + f^rY — + ^ [93.75, 125] GeV is applied to emulate the 
fit to the distribution performed in the ATLAS analysis. 

^We implement the relevant quality and isolation criteria at the level of the Delphes simulation. 
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Fig. [TT] displays the m« (upper panels) and 0« (lower panels) distributions for Xo+_o- (left 
and centre panels) and X2 (right panels) after implementing in Delphes the ATLAS analysis 
cuts described above. We see that the effect of cuts is dramatic, reshaping the distributions 
in the X2 case so that it resembles the Xq hypothesis. This is not only a consequence of 
the A(j)u cut. We have verified that one could loosen or even remove the A0« cut: its effect 
on the background rejection is very mild, and dropping it would not help to maintain the 
distinctive kinematic features of X2 — )■ WW* — )■ i'^uii'ui decay. The initial cut plays a 
key role in reducing features due to the anti-parallel preference of the lepton momenta in 



Fig. 11, due in turn to the strong correlations between the Acpu, m« and cuts. 
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Figure 11: Results from simulations of the Xo+_o-,2+ ~^ WW* — )■ C^vil'vi signals, using 
PYTHIA and Delphes and implementing the ATLAS cuts described in the text. The upper 
(lower) panels display themu (4>u) distributions for Xq+ fl- (left and centre) and X2+ (right). 



In Fig. 12, we show the distributions found after implementing all the ATLAS cuts 
described above, under the Xq+^q- and graviton-like X2+ hypotheses. They exhibit somewhat 
different behaviours within the selected kinematic range. The peaking of the X2+ histogram 
at slightly lower rriT than those for Xq+^q- is a relic of the differences in the kinematic distri- 
butions before the cuts. Since the neutrino and and antineutrino are emitted antiparallel in 
X2+ case, the invariant mass has a slight tendency to be larger than in the Xq+^q- cases, 
implying that rriT tends to fall further below the X mass of ~ 125 GeV. 
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Figure 12: The transverse mass, ttlt, distribution after application of all the ATLAS cuts 
described in the text, under the three hypotheses. 

We have also simulated the corresponding CMS search for Xo,2 — > WW* — > l^Pil'vi. 
The CMS cuts are very similar to those applied by ATLAS, except that CMS requires > 
45 GeV, mu e [12,45] GeV, A^ee < 1-6 and tjit e [80, 125] GeV. The resulting histograms 
of mu and (pu are very similar to those for ATLAS shown in Fig. 11 , so we do not show the 
CMS equivalents. 

We display in Table [T] the ATLAS and CMS cut flows for the different X jp hypotheses. 
These results are based on simulations of 10,000 0^ events, 20,000 0~ events and 30,000 
2"*" events, so the statistical errors are negligible. Although there are differences between 
the numbers of events surviving different stages in the ATLAS and CMS selection, the end 
results after applying all the cuts are very similar. Specifically, we find that the efficiencies 
for the O"*" hypothesis are larger than those for the 2"*" hypothesis by factors of 1.86 (1.94) 
for the ATLAS (CMS) event selections, with the efficiencies for a pseudoscalar Xq- being 
about 10% lower than that for Xq+ in both cases. 



4.2 Data analysis under different hypotheses 

We now discuss how this efficiency difference in the experimental selection cuts may in prin- 
ciple be used to help discriminate between the scalar and graviton-like spin two hypotheses, 
analyze the sensitivity offered by the current data and estimate the likely sensitivity of the 
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Lilts 


bCfVirVl (J^ 


Psoudoscalar 


LiaA lt()ll-ilkC> 


pi^ > 30 (45) GeV 


86 (63) 


83 (58) 


70 (43) 


m^^ < 50 (45) GeV 


66 (49) 


62 (44) 


40 (26) 


A(!>u < 1-8 (1.6) 


63 (49) 


60 (44) 


38 (26) 


rriT e [93.75, 125] ([80, 125]) GeV 


44 (40) 


41 (36) 


22 (21) 



Table 1: The outflow evolutions for the ATLAS (CMS) cuts under the three hypotheses. 
The numbers shown are the cumulative efficiencies to pass the cuts (in percent). 

full 2012 data set. 

We parametrize the rescaling of Xo+ particle couplings to the W, Z gauge bosons relative 
to the Standard Model values by aw and az, respectively, and infer Xwz = dw/^z from the 
measured ratio of the signals in the WW* and ZZ* channels, which is given by 



R 



wz 



Nw BRsm{Xq+ WW* 



az 



(17) 



BRsm{Xo+ ^ ZZ*) 

where ew,z are the efficiencies for the Xo+ — )■ WW*, ZZ* experimental selections. Likewise, 
rescaling by aw2 and a^^ the W and Z couplings of an X2+ particle relative to reference 
values with custodial symmetry, one has in an obvious notation 



R 



BR2+iX2+ WW*) [aw, 



~ Nz ~ BR2+{X2+ ^ ZZ*) \az2J \ez2J' ^ ^ 

which can be used to infer Xwz2 = Ciw2/^Z2 i^i the same way. Since we use X — >■ ZZ* — )■ 
event selections that use only the individual momenta (specifically, we do not use the 
CMS MELA analysis), we may assume that e^j/e^ = 1. The value of \wz2 inferred from 
the data therefore differs from that of Xwz by the following factor 

/ BRsMiXo+ ^ WW*) BR2+iX2 



^WZ2 — ^wz X 




Co+ ^ ZZ*) ' BR2+{X2+ ZZ* 



(19) 



where the value ew/^W2 — 1-9 was calculated in the previous Section. 

The ratio of the ratio of Xo+_2+ — > WW* and ZZ* branching ratios is not simply unity, 
because of the non-trivial dependences of the partial decay widths T{Xq+^2+ — > VV*) on 
the masses of the vector bosons V. We have used MadgraphS vl.4 and vl.5 to calculate 
the decay widths, and have checked our results against the Standard Model predictions for 
H — > WW* and ZZ*, and also tested them in the limiting cases of heavy graviton and Higgs, 
when the vector bosons are produced on-shell. In the case of the physical X mass, we find 

BR2+{X2+ WW*) .BRsm{Xo+ ^ WW* 



BR2+{X2+ ^ ZZ*) ' BRsm{Xo+ ZZ* 



= 1.3, 



(20) 
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ATLAS 7+8 WW.ZZ ATLAS 7+8 WW.ZZ 




Figure 13: Left: The {aw,0'z) plane displaying the 68 and 95 % CL contours (broken black 
and solid blue, respectively) obtained from a fit to the ATLAS 7 and 8 TeV data, with rays 
corresponding to contours of aw/az: the Standard Model point aw = az = I is indicated 
by a black star. Right: The Ax^ function relative to the best-fit value in our ATLAS fit, 
marginalized over the magnitudes of aw,z, with dashed horizontal lines at Ax^ = 1,4. 



and hence 

^WZ2 = ^wz X 1-2, (21) 

when aw and az, and hence Xwz and Xwz2, are extracted using only the X — )■ WW* and 
ZZ* inclusive search channels. 



In Fig. 13 we display in the left panel the {aw, az) plane with the two-dimensional CL 
contours that we find in our analysis of the ATLAS 7 and 8 TeV data combined |3j . Also 
shown to guide the eye, in this and subsequent similar plots, are rays corresponding to various 
values of the ratio aw/az- The ray aw/az = 1 passes through the Standard Model point 
aw = az = ^, which is indicated by a black star. This point also lies just within the 68% 
CL contour, shown as a broken black line (the 95% CL contour is a solid blue line). The 
right panel of Fig. 13 displays the Ax^ function relative to the best-fit value in our analysis 
of the ATLAS 7 and 8 TeV data, marginalized over the magnitudes of aw,z- The combined 
ATLAS data do not exhibit any strong preference between the = 0+ and 2"*" hypotheses, 
which correspond to aw /az = ^ and 1 / \/2, respectively. 

Fig. 14 displays a similar pair of panels for our analysis of the available CMS 7 and 8 TeV 
data. We note here that we do not use the final CMS result for X — )■ WW* signal, which 
include an MELA selection that we do not model. Instead, we use the expected signal, 
background and observed event numbers shown in Table 3 of [2j, which correspond directly 



We cannot use for this purpose the value of Xwz quoted in 24 , because this incorporates information 
from a combination of channels including X — ^ 77, which supplies information on aw that does not apply 
to the spin-two case. 
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Figure 14: As in Fig. 13_, hut based on our fit to the CMS 7 and 8 TeV data, which does 
not include the final MELA selection made in 



to the CMS event selection and efficiency found in our simulation above. As in the ATLAS 
case above, the CMS data also do not exhibit a strong preference between the = 0+ and 
2"*" hypotheses. 



Fig. 15 displays a similar pair of panels for our combined analysis of the available AT- 
LAS and CMS 7 and 8 TeV data. The data sets in combination provide no significant 
discrimination between the = 0"^ and 2"*" hypotheses. Our combined analysis yields 



A 



wz 



az 



u.yo_o.2o ' 



Using (21), the corresponding value in the X2+ case is 



A 



WZ2 



1 19+0.29 

-'-•-'-^-0.24 1 



(22) 



(23) 



and we see that both results are compatible with unity. 

This result is based on ~ 5/fb of data at each of 7 and 8 TeV analyzed by each of 
ATLAS and CMS. At the time of writing, each experiment has now recorded ~ 15/fb of 
data at 8 TeV, corresponding approximately to a doubling of the statistics and potentially 
to a reduction in the uncertainty in a^/jaz by a factor ~ 1.4. It is anticipated that each 
experiment might record ~ 30/fb of data by the end of 2012, corresponding to reductions 
in the statistical errors in (22, 23) by factors ~ v^- If the central value remained the same. 



(23) would become \wz2 = ^■^'^-o'll^ whereas if the central value were to correspond to 



^wz = I5 one would have Xwz2 = l-21_ai4- We conclude that the best one could reasonably 
hope for with the 2012 data would be a deviation from the custodial symmetry prediction 
Xwz2 = 1 of about one a. 
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Figure 15: As in Fig. hut based on our combined fit to the ATLAS and CMS 7 and 
8 TeV data. 



5 Overview and Prospects 

We have explored in this paper two strategies to help determine the spin of the new particle 
X discovered recently by the ATLAS and CMS Collaborations. One of these exploits the 
angular distribution of the final-state photons in (7(7 — )• X -> 77 decay, and the other 
exploits the angular distributions and correlations in X — )■ WW* — )■ t^t~vv decay, which 
are in principle quite different for different spin assignments = 0~^, 0^ and 2+ for the X 
particle. 

We have shown how the 2012 LHC data could be used to study the angular distribu- 
tion in (7(7 — > X — 77 decay and provide potentially significant discrimination between the 
spin- zero and spin- two hypotheses. A simple angular asymmetry measurement gives a dis- 
crimination power approaching that possible with a full LLR analysis. We include the effects 
of backgrounds in samples with both high and low S/B = 0.42 and 0.19, corresponding in 
spirit to CMS event categories. We find that the present data should already provide some 
discrimination between the = 0^ and 2+ hypotheses, and that an analysis of the full 2012 
data set could provide a separation significance of ~ 3 cr if a conservative symmetric inter- 
pretation of the LLR statistic is used, rising to above 6 cr if a less conservative asymmetric 
interpretation is used. 

We have analyzed the sensitivities of the published results of the ATLAS and CMS 
searches for H — )■ WW* to the spin-parity of the X particle. Simulating these searches using 
PYTHIA and Delphes, we have shown that the ATLAS and CMS experimental selections 
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suppress the kinematic differences between 0"^ and 2+ decays. Therefore, an analysis based 
on kinematic shapes would be rather inconclusive based in the published cuts. 

One could hope for retaining some of the kinematic differences by changing the cuts, but 
those changes are hmited by the need to suppress the large Standard Model WW background. 
A more hopeful strategy, and the one we developed in this paper, is to use of the approximate 
custodial symmetry in both spin-zero and -two hypothesis, to relate the WW* and ZZ* 
channels. In WW*, the efficiencies of the searches differ by a factor ~ 1.9 for Xo+ and X2+. 
On the other hand, we find that the ratio of X2+ — > WW* and ZZ* branching ratios is 
~ 1.3 larger than the corresponding ratio of branching ratios in the 0+ case. The current 
measurements by ATLAS and CMS of the ratio of experimental rates for X — > WW* — > 
t^l'vV and X — >■ ZZ* — >■ M'^ are compatible with custodial symmetry in both the 0+ and 
2+ cases, and we do not expect that this will improve significantly with the full 2012 set of 
LHC data. 

Many other strategies to discriminate between different spin-parity hypotheses for the X 
particle have been proposed, including kinematic correlations in X — > ZZ* — > 4^=*= decay and 
the threshold behaviour of associated W/Z-\-X production, as discussed in the Introduction. 
We have every confidence that the spin-parity of the X particle will be pinned down within 
a few months, in good time for the nominations for the 2013 Nobel Physics Prize. 
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